Gas dynamics in high-luminosity polarized ^He targets using diffusion and convection 
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The dynamics of the movement of gas is discussed for two-chambered polarized ^He target cells 
of the sort that have been used successfully for many electron scattering experiments. A detailed 
analysis is presented showing that diffusion is a limiting factor in target performance, particularly 
as these targets are run at increasingly high luminosities. Measurements are presented on a new 
prototype polarized ''He target cell in which the movement of gas is due largely to convection instead 
of diffusion. NMR tagging techniques have been used to visualize the gas flow, showing velocities 
along a cylindrically-shaped target of between 5 — 80 cm/min. The new target design addresses one 
of the principle obstacles to running polarized ^He targets at substantially higher luminosities while 
simultaneously providing new flexibility in target geometry. 

PACS numbers: 29.25.Pj, 13.60.Hb, 13.60.Fz, 25.30.Bf 
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I. INTRODUCTION 

Nuclear-polarized ^He has proven to be useful in a 
number of different areas of research. In electron scatter- 
ing, polarized ^He provides a means for studying spin- 
dependent interactions involving neutrons. This is be- 
cause, to first approximation, a "^He nucleus is comprised 
of two protons whose spins are paired, and a single neu- 
tron that accounts for most of the nuclear spin [1^ . An 
important early example of the use of polarized ^He in 
electron scattering came during an experiment to mea- 
sure the internal spin structure of the neutron at the 
Stanford Linear Accelerator Center (SLAC), E142 [5]. 
Polarized '^He has also been used to measure the electric 
form factor of the neutron G^, including a recent exper- 
iment at Jefferson Laboratory (JLab) in Newport News, 
VA [3,. Important applications of polarized "^He have also 
included its use as a neutron polarizer [3], and, together 
with polarized "'^^^Xe, as a source of signal for magnetic 
resonance imaging |5l E] . 

There are two predominant techniques by which high 
levels of nuclear polarization are produced in ^He. In one 
technique, often known as metastability exchange opti- 
cal pumping (MEOP), metastable states of ^He are opti- 
cally pumped directly, and subsequently transfer their 
polarization to other ground-state '^He nuclei during 
metastability- exchange collisions [3 [S] . In a second tech- 
nique, known as spin-exchange optical pumping (SEOP), 
a vapor of alkali-metal atoms is optically pumped and 
subsequently transfers its polarization to '^He nuclei via a 
hyperfine interaction during spin-exchange collisions [S]- 
[TT] . An important difference between the two techniques 
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is that MEOP is performed at pressures that are quite 
low, around a few Torr, whereas SEOP is often done at 
pressures as high as roughly ten atmospheres. When a 
high-density is required, a target based on MEOP in- 
evitably involves a compressor of some sort. In contrast, 
high-density targets based on SEOP typically utilize a 
sealed glass cell with no moving parts and hence have an 
advantage from the perspective of simplicity. In consid- 
ering relative merits, however, one also needs to consider 
the speed with which the gas is polarized. Here targets 
based on MEOP have done quite well, with a recent tar- 
get at the Mainz Microtron reporting a polarization rate 
of 2 bar-liters per hour|12J. In short, both techniques 
for polarizing ^He have been quite successful and have 
complimentary advantages. 

For electron-scattering experiments in which "^He is po- 
larized using SEOP, the targets typically utilize a sealed 
glass cell with two distinct chambers: a "pumping cham- 
ber" in which the gas is polarized, and a "target cham- 
ber" through which the electron beam passes (see Fig. II]) . 
This design ensures that ionization due to the electron 
beam does not adversely affect the optical pumping pro- 
cess, as well as providing in the pumping chamber a ge- 
ometry that lends itself well to illumination with lasers. 
The two chambers are connected by a "transfer tube", 
and gas that is polarized in the pumping chamber mi- 
grates into the target chamber largely by diffusion. 

Several conditions need to be met in targets with de- 
signs such as that shown in Fig. 1 in order to maintain 
high polarization. First, the rate at which the "^He is 
polarized must be relatively fast compared to the rate 
at which the '^He is depolarized. While this is true in all 
work involving SEOP, it is of particular significance when 
considering depolarization of the '^He due to an electron 
beam. There are also issues having to do with polariza- 
tion dynamics that are unique to the two-chamber de- 
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FIG. 1: Shown is the geometry of a two-chambered glass ceU 
used for polarized ^He targets. The dimensions shown are 
typical of those used in a recent JLab experiment to measure 
the electric form factor of the neutron |3] • 



sign. As has been pointed out by Chupp et al.[13\, it 
is important that diffusion between the two chambers is 
rapid compared to other '^He-related time constants in 
the system. As we will show in detail, a matter of special 
importance is maximizing the ratio of the diffusion rate 
to the relaxation rate that is specific to the target cham- 
ber. Put differently, it is critical that the polarized gas in 
the target chamber is replenished much faster than it is 
depleted through depolarization mechanisms such as, for 
example, ionization from the electron beam. A failure to 
replenish the polarized gas quickly enough will result in a 
lower polarization in the target chamber than is the case 
in the pumping chamber, a condition we refer to herein 
as a polarization gradient. 

Until relatively recently, polarization gradients in two- 
chambered cells have been at most a few percent rela- 
tive. Advances in SEOP, however, have made it possible 
to run polarized '^Hc targets at increasingly high lumi- 
nosities. During SLAC E142, where luminosities were in 
the range of 0.42 — 1.70 x 10'^^ cm~^s~^, the polarization 
gradients were on the order of 1%. During recent experi- 
ments at JLab, however, with luminosities in the range of 
0.6 — 1.0 X 10'^^ cm~^s~"'^, polarization gradients were as 
high as nearly 9% (relative) . For future experiments with 
luminosities in the range of lO'^^cm^^s^^, polarization 
gradients well in excess of 10% are likely if no changes 
are made to the basic target-cell design. Polarization 
gradients can also be difficult to quantify accurately, an 
issue that can lead to uncertainties in polarimetry. 

The fact that polarized ■^He targets are being run at 
increasingly high luminosities is due largely to advances 
in SEOP. One example is the use of hybrid mixtures of 
alkali metals (typically potassium and rubidium) instead 
of a single alkali metal (typically rubidium) [TH [TS] . This 
technique, often known as alkali-hybrid SEOP, greatly 
improves the efficiency with which the '^He is polar- 
ized. Another important advance has followed from 
the availability of commercial spectrally-narrowed high- 



power diode-laser arrays. These new lasers result in 
significantly higher optical pumping rates for a given 
amount of light. Collectively, these advances have made 
it possible to significantly increase the rate at which the 
■^He is polarized. It has thus become possible to achieve 
higher polarizations even when using higher beam cur- 
rents and thicker targets. Higher currents, however, make 
it necessary to replenish the polarized gas in the target 
chamber more quickly. If the advances in SEOP are to be 
fully exploited, it is essential that the designs of polarized 
■^He targets evolve. 

We report here on the successful implementation of a 
new design for polarized '^He target cells based on SEOP. 
The design incorporates the ability to circulate the gas 
between the pumping chamber and the target chamber 
using convection instead of diffusion, an idea discussed 
by Wojtsekhowski in Rcf. [16 in anticipation of both 
advances in SEOP as well as the need for higher lumi- 
nosities. The convection is achieved by maintaining a 
temperature differential between different parts of the 
target cell, and does not involve pumps or other moving 
parts. We have shown that the velocity of the gas mov- 
ing through the target chamber can be varied between 
5-80 cm/min in a simple controllable manner. The ad- 
vent of a means to circulate a polarized noble gas in a 
sealed vessel without the use of pumps has great potential 
for high-luminosity polarized '^He targets. The simplic- 
ity of the approach has advantages from the perspective 
of reliability, and fast transport of gas between the two 
chambers makes it possible to greatly increase the elec- 
tron beam current without causing a polarization gradi- 
ent. Convection-based target cells also open the possi- 
bility of physically separating the pumping chamber and 
the target chamber by much larger distances than was 
previously possible, something that offers several practi- 
cal advantages. 



II. POLARIZATION DYNAMICS AND 
GRADIENTS IN TWO-CHAMBERED CELLS 

A. The Single-Chambered Cell 

Before considering the formalism for target cells with 
two chambers, we begin by considering the simpler ex- 
ample of a single-chambered cell, where the equation de- 
scribing the time evolution of the polarization is given 
by: 



Ate = 7so-Pa - (7sc + r)PHc, 



(1) 



where Phc is the '^He polarization. Pa is the polarization 
of the alkali vapor, 7se is the rate at which the '^He is po- 
larized due to spin exchange, and F is the spin-relaxation 
rate of the '^He due to all other processes. The solution 
to Eqn[l]is given by 
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where Pq is the ^He polarization at i = 0. It has been 
shown by Babcock et al. that one of the components 
of r is a relaxation rate that empirically appears to be 
proportional to 7se |17| . One can accordingly write that 
r = r' + 7sG^, where X is a proportionality constant. 
We can thus write the saturation polarization associated 
with Ecyi. [2] as 
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where we note that the denominator of Eqn. [3]is also the 
rate that characterizes the buildup of polarization in [2] 
rewritten in terms of T' and X. The existence of a relax- 
ation mechanism proportional to 7se is unfortunate, as it 
implies that, even with a 100% polarized alkali vapor, in 
the limit that 7so — >■ 00, Pho -^ ^/i^+X). While not well 
understood, the existence of this additional relaxation is 
now well established, and is important for understanding 
the overall polarization achieved. 

It is straightforward to project the performance of a 
target at an arbitrary range of beam currents if we know 
how it performed at a particular beam current /q (true 
even if Iq — 0). Let us assume that we have data describ- 
ing polarization as a function of time while polarizing a 
target from an initial polarization of zero. A plot of data 
of this sort is something we refer to herein as a "spin-up 
curve", and can be fit to yield the "spin-up rate" that 
appears in the denominator of Eqn. a\ (as well as the 
exponent in Eqn. p] 7sc + T)- In keeping with the nota- 
tion introduced above with respect to Eqn. [3) we define 
a quantity 7°^ that characterizes the spin-up rate at a 
particular current Iq: 
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Let us further define -Phc(-^o) as the equilibrium polariza- 
tion associated with that spin-up. Using Eqns. |3] and |4j 
we find 
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The quantities P^^ilo) and 7"^ can be determined by 
fitting data from a spin up, and the quantity r'(/) — 
Twin + Tdipoio + Lbcam, whcrc the three terms are spin- 
relaxation rates due to wall collisions (at room tempera- 
ture), dipole interactions due to ■^He-'^He collisions, and 
ionization of the electron beam respectively. The sum 
of the first two terms is essentially the target's room- 
temperature spin-relaxation rate in the absence of beam, 
a quantity that is quite easy to measure. We note, how- 
ever, that Fdipoio will be slightly different under operating 
conditions than it is at room temperature both because of 
heating the pumping chamber and the related differences 
in densities. This correction is easily calculated using the 
calculation of Newbury et al. from Ref. [18] . The contri- 
bution Fboam is also easily computed [HI |2Q] . We note 
that calculations of Fbeam have shown good agreement 
with experiment at a level of roughly 10% or better [2T] . 



We can thus use Eqn. [5] to project the performance of a 
particular target at an arbitrary beam current. 

It is instructive to investigate how different targets that 
have been used during past experiments would fare at 
significantly increased current, for example, at 100 /iA. 
Here we ignore the effect of polarization gradients. The 
first time liter-type quantities of polarized ^He were used 
in a target was the aforementioned experiment at SLAG, 
E142 [2]. In the presence of 3.3 /iA of electron beam, 
the '^He polarization averaged about 33%, with values 
for (7°u)~"^ of about 15-20 hours. If this same target 
were instead exposed to 100 /iA of electron beam, Eqn. [5] 
suggests that the '^He polarization would drop to just 
over 10%. In contrast, during the more recent experi- 
ment at JLab that measured G^ |3], the cell-averaged 
polarizations were around 50% with 8 /iA of beam and 
(7su)~^ was in the range of 5 — 6hrs. Here Eqn. ^sug- 
gests that at 100 /iA the resulting "^He polarization would 
be around 37-38%. The improved projected performance 
is due both to the shorter values for (7°u)~^j as well as 
the fact that the target cells were significantly larger, 
making them more resistant to depolarization from the 
electron beam. We note, however, that even though the 
cell-averaged polarization would be fairly reasonable, the 
polarization that one would have in the target chamber 
would be much lower because of polarization gradients. 



B. Time Evolution in a Double-Chambered Cell 

For a full description of a double-chambered cell, the 
polarization build up must be described by the coupled 
differential equations first described in Ref. [13| : 
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^tc — ^Ftc^tc + "tc(^pc ~ ^tc) 
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where Ppc (Ptc) is the ^He polarization in the pumping 
(target) chamber, 7so is the spin-exchange rate in the 
pumping chamber, and Fpc and Ftc are the ^He spin- 
relaxation rates due to all other processes in the pumping 
and target chambers respectively. The transfer rate dtc 
(dpc) is the probability per unit time per nucleus that a 
nucleus will exit the target (pumping) chamber and enter 
the pumping (target) chamber. We note that we do not 
include the transfer tube as a separate volume in this 
analysis. The transfer rates are related by 



/pc'^pc — /tc^tt 



(8) 



where fpdftc) is the fraction of atoms in the pumping 
(target) chamber, and /pc + ftc = 1- For the dynamic 
studies reported in Ref. T3], the authors were able to 
neglect terms involving 7sc and F relative to terms in- 
volving dpc and dtc. For the discussion here, however, we 
must retain these terms, requiring an analysis essentially 
identical to that considered by Jones et al. [22] and later 
by Kominis [23]. We refer the reader to those two refer- 
ences for details. We find that the solutions to Eqns. [6] 



and [7] are given by 

Ppc(i) = Cpee-i^'* + (Pp", - Pp- - Cpc)e-r=* + P- (9) 
and 



the same form as the time constant that appears in Eqn.[2] 
(7sc + r). It can be shown that 
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where Pp^. and P^° are the initial polarizations in the 
pumping and target chambers respectively, 
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We liave chosen to write Eqn. [TT] in the form shown to 
emphasize that 
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where (7sc) is the spin-exchange rate averaged through- 
out the double-chambered cell ((7so) = /pc7so, since the 
spin-exchange rate is 7se in the pumping chamber and is 
zero in the target chamber), and (F) = /pcFpc -I- /tcFtc 
is the spin-relaxation rate averaged throughout the cell. 
Eqn. [13] has the same form as the saturation polarization 
of Eqn. [2J as we would expect in the limit of infinitely 
fast transfer. 

The coefficients Cpc and Ctc are given by 
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where a = -{%c + Tpc + dpc), b = dpc, c = dtc and 
d = —(Etc + dtc)- We note that in the fast-transfer limit, 
the quantities Cpc and Ctc are given by 



and 



Cpc — /tc(Ppc ~ ^tc) 



Ctc — /pc(Ptc - ^pc) 
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The constants Eg and Ef represent slow and fast rates 
respectively that govern the time evolution of the polar- 
ization. It is useful to write E, in the form 



(7sc> + (r> - ^E 



(18) 



where the quantity ST is generally small and goes to zero 
in the limit of infinitely fast transfer between the two 
chambers. We can see that the limiting form of Eg has 



7sc + Tpc - Et, 

"pc + "tc 



(20) 



Eor most of the situations we would normally consider, 
the quantity u is fairly small. This is due to two things. 
Eirst, the spin-exchange rate 7se is typically slow com- 
pared to the sum of the two transfer rates dpc and dtc, 
and second, both Epc and Etc rnust be relatively small 
compared to 7so, or the polarization of the target would 
not be high. It is thus reasonable to expand 6T in a 
Taylor series in u: 

(^r sa fpcftc{dpc + dtc)u^ + higher order terms. (21) 

Lastly, we consider Ef which can be written as 

Ef = (dpc + dtc) + (7sc - (7sc)) + (Fpc + Etc - (E)) + <5E . 

(22) 
In the fast-transfer limit, Ef — > 00; under these condi- 
tions, Eqns. [9] and 10 reduce to the form of Eqn. [2] 

Data illustrating the time evolution of polarization 
(what we referred to earlier as a spin-up) are shown in 
Eig. [2] for both the pumping and target chambers of a 
double-chambered cell we refer to herein as "Brady" . The 
polarization was measured every three minutes using the 
NMR technique of adiabatic fast passage (AEF) [21]. We 
note that under normal operating conditions, NMR mea- 
surements would only be made once every few hours, in 
part because each measurement results in a small loss 
(< 1%) of polarization. The frequent measurements 
shown in Eig. [2] strongly limit the saturation polarization 
because of accumulating losses. Also shown in Eig. [2] 
but obscured beneath the many data points, is a fit to 
the data using double-exponential functions of the form 
given in Eqns. |9] and 10 The fit clearly describes the 
data quite well. 

Einally, we note that in the context of the types of cells 
that have been used in electron-scattering experiments 
(see Fig. II]) , the mechanism behind the transfer rates dtc 
and dpc is overwelmingly diffusion. 



C. Initial Polarization Evolution 

Some of the parameters discussed earlier can be readily 
determined by studying spin-up curves of the sort shown 
in Eig. [2| To extract values for the transfer rates dp^. 
and dtc, it is particularly valuable to examine the spin- 
up curves for the initial time period during which the 
polarization is growing. Eor small values of the time t, it 
is readily apparent from Fig. [2]that the nature of the time 
evolution in the two chambers is quite different. Under 
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FIG. 2: The ^He polarization is shown as a function of time for 
both the pumping chamber (upper curve) and target chamber 
(lower curve) of the target cell "Brady" . In this case the lasers 
were turned on immediately before data taking to ensure an 
initial polarization of zero. Also shown are fits to Eqns.joland 
|10| We refer to curves of this sort as spin-up curves. AFP 
measurements were made rapidly (every 3 minutes). 



the assumption that the time t < l/Ff (in this case, 
l/Ff sa 0.75hrs), we can expand Eqns. [9] and 10 
Taylor series. ~ . . - - - — 



To second order, for the case of P^ 



0, this expansion simphfies to: 
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and 
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In Fig. [3j we show only the first 24 minutes of the data 
shown in Fig. [2] It can be seen that the initial shape of 
the spin-up curve appears to be linear in the pumping 
chamber and quadratic in the target chamber, in agree- 
ment with expectations from Eqns. [23] and [24j 

To empirically determine dtc, we note first that the 
slope of the nearly-linear polarization buildup in the 
pumping chamber is equal to the product Pa7so- With a 
fit to this slope, along with a fit to the coefficient charac- 
terizing the quadratic polarization buildup in the target 
chamber, we can extract a value for the transfer rate 
dtc = (0.72±0.10)hrs~^. As will be discussed in the next 
subsection, a value for dtc can also be computed from first 
principles given the dimensions of the cell. The compar- 
ison of empirically determined and calculated values for 
dtc provides insight into our understanding of the diffu- 
sion processes taking place in our cells. 



D. Transfer Rates Under Diffusion 

Using gas kinetic theory, the dimensions of the target 
cell, the density of '^He and other gases when the cell was 



filled, the operating temperatures of the pumping and 
target chambers and the assumption that the tempera- 
ture gradient along the transfer tube is linear, it possible 
to compute dpc and dtc from first principles. We begin 
similarly to the discussion in Ref. [T31 with an equation 
describing the net polarization flux Jtt through the trans- 
fer tube: 
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dz 



(25) 



where n{z), D{z) and P{z) are the '^He number density, 
self-diffusion constant and polarization respectively, all 
shown as a function of position z along the transfer tube. 
As discussed by Romalis ^ and later by Zheng \2& , us- 
ing data on the self-diffusion constant of '^He by Kestin 
et al. [27], the self diffusion constant for ■^He is well ap- 
proximated by the expression 
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= 353.14 K, m = 1.705 and 
2.687 X lO^^cm-^). The 



where Dq = 2.789 cmVs, To 

no = 0.7733 amg (1 amg = 

assumption of a linear temperature gradient along the 

transfer tube between the pumping and target chambers, 

as was assumed in Refs. [25, 26J, allows us to express 

T{z) and hence n{z) explicitly. With this assumption. 



and substituting Eqn. J26] into Eqn. 25 we can solve for 
Jtt to find 
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Jtt is the total rate of polarization transfer per unit area, 
whereas we want the rate per atom. Multiplying by the 
transfer tube cross sectional area Att, dividing by the 
number of particles in each chamber, and finally dividing 
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FIG. 3: Early time behavior of the Brady spin-up shown in 
Fig.i 



by (Ppc - -Ptc), we have 



1. Quantifying polarization gradients 
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Eqn. [28 
0.667cm 
11.5 amg 



Of specific interest here is the value for dtc implied by 
for the target cell Brady. For Brady, Att — 
, Ltt — 9.07cm, Vtc = 74.6 cm'^, and rite = 
Using temperatures that correspond to the 
tests iUustrated in Figs.[2]and|3J we find dtc = 0.72 hrs~^, 
a value that agrees fortuitously with the value found by 
fitting the polarization buildup curves. 



E. Polarization Gradients 

An issue of considerable practical importance for po- 
larized '^He targets is the polarization gradient between 
the pumping and target chambers. Dividing both sides 
of Eqn. 
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by P^, we find 
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It is also useful to define the quantity A, the amount 
by which the polarization in the target chamber is lower 
than that of the pumping chamber. Here we express this 
difference as a fraction of the pumping chamber polar- 
ization: 
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We can see that A approaches for cells in which Ftc is 
quite slow and dtc is quite fast. The fact that it is Ftc 
and not the cell-averaged relaxation rate (F) that appears 
in this expression is quite important. When a target is 
subjected to high electron-beam currents, the overall cell- 
averaged relaxation rate (F) may not be strongly affected 
even though the local relaxation rate in the target cham- 
ber Ftc is. As stated earlier, what is important is that 
the polarized gas in the target chamber is replenished at 
a rate that is much faster than the rate at which the gas 
is depolarized in the target chamber. 

Looking at specific examples, we find that the targets 
used during E142 at SLAG had polarization gradients 
of just over 1% with no beam current, and would hy- 
pothetically have had polarization gradients of around 
13% at 100 /xA. For the targets used to measure G^ at 
JLab, the gradients were worse, a little over 6% with no 
beam current, and we project they would be about 21% 
at 100 /iA. The larger gradients in the case of the G^ cells 
were mostly due to faster intrinsic cell-relaxation rates, 
but also to having longer transfer tubes, something that 
is difficult to avoid with a cell that has a larger overall 
size. Shortly we will return to the question of what can 
be done to better minimize polarization gradients when 
using a single transfer tube. 



Regardless of the magnitude of the polarization gra- 
dient, it is also difficult to quantify accurately because 
of uncertainty in the quantity Ftc. When characteriz- 
ing a target cell, the quantity that is most straightfor- 
ward to measure is the cell-averaged room-temperature 
spin- relaxation rate (F), which is typically due to three 
primary contributions: 



(F> = F" 



(31) 



where F™ is the spin relaxation rate due to wall collisions, 
F^ is spin relaxation due to dipolar interactions during 
•^He — ^He collisions and F* is the spin relaxation due 
to the electron beam, which can be taken to be zero if 
the cell's spin-relaxation rate is measured in the absence 
of an electron beam. Here we ignore relaxation due to 
magnetic field inhomogeneities which with some care can 
be made quite small. The wall relaxation rate will be 
the sum of the wall relaxation rates in the target and 
pumping chambers respectively, weighted by the fraction 
of '^He atoms that are in each chamber: 



/tcrtc 



Jpc^ pc 



(32) 



For the purposes of this discussion, it is convenient to 
introduce a parameter R, representing the ratio of F^ to 



F^^, so that 



J^ tc — -'^ J^ pc 



(33) 



(30) From Eqns 



31 



33 



taking F'' = 0, we find an expression 



for F^ 



^((r)~F'^) 

/tc R + /pc 



(34) 



Unfortunately, we have no direct measurement of R, and 
wall-relaxation rates are notoriously variable. One plau- 
sible assumption is that F^ — Fp^, in which case R = 1, 

and Eqn. [sllsimphfies to Ft^c = Tpc = (r) - r"*- This is. 



in fact, the assumption that has been made in polarized 
•^He experiments (those using the basic design shown in 
Fig. IT]) prior the aforementioned G^ experiment. Sev- 
eral authors have shown, however, that for uniform wall 
relaxation per unit area, overall wall relaxation is propor- 
tional to the surface-to-volume (S/V) ratio of the vessel 
containing the gas [211 [H]. This dependence was recently 
exploited by Anger et al. who successfully constructed 
storage cells for polarized ^^^Xe with unusually long, per- 
haps unprecedented, relaxation times [30^ . If wall relax- 
ation were uniform throughout a target, we would expect 
R to be the ratio of the S/V ratios of the pumping and 
target chambers respectively, a quantity we will refer to 
here as i?max- A conservative approach, then, might be to 
take R — {1 -|-i?niax)/2, with an error on A that includes 
the full range of 1 < ii < 7?max- 

The uncertainties in Ftc, and consequently A, can in 
some situations translate into a systematic uncertainty 



in polarimetry. One of the best techniques for deter- 
mining the absolute polarization of "^He is the method of 
measuring shifts in the electron paramagnetic resonance 
frequencies of the alkali-metal atoms due to the effective 
magnetic field caused by the polarized '^He [31 . This 
can only be performed in the pumping chamber (where 
significant alkali-metal vapor is present), despite the fact 
that the quantity of interest in an electron scattering ex- 
periment is the polarization in the target chamber. Thus, 
it is often the case that NMR measurements are made di- 
rectly on the target chamber, but are calibrated against 
frequency shift measurements in the pumping chamber. 
When this is done, the calibration requires a knowledge 
of A. For the case of the targets used in the G^ exper- 
iment discussed earlier, the uncertainty in A, following 
the prescription described earlier, translated into a 1-2% 
(relative) uncertainty in polarimetery. While not catas- 
trophic, such systematic uncertainties would be nice to 
avoid. 



as large as is practically realizable). We further assume 
the transfer tube has the same length as before, and that 
the pumping chamber is larger (around 12.5 cm) so as to 
satisfy the criterion of having the target contain 7 STP 
liters of gas. Such a cell would have a polarization gradi- 
ent of around 12% in 60 /xA of beam. If the cell-averaged 
polarization were around 70% with no beam (more on 
this in section IV), and 62% in beam, the polarization in 
the target chamber would be about 55%. In addition to 
the reduction of the in-beam polarization, the gradients 
would also have the potential of introducing uncertainties 
in polarimetry at the level of 2-4%. In short, this hypo- 
thetical design falls well short of optimized performance. 
And while we have assumed here the same transfer-tube 
length as is shown in Fig. 1, as will be discussed more in 
section IV, there are compelling reasons to increase the 
distance between the pumping and target chambers. This 
would naturally require a longer transfer tube, something 
that would further aggravate the problem of polarization 
gradients. 



2. Limitations from, using a single transfer tube 

Two points emerge regarding polarization gradients: 
1) they limit target polarization at high beam currents, 
and 2) they are somewhat uncertain in their size, which in 
some circumstances, results in systematic uncertainties in 
polarimetry. If we want to make polarization gradients 
smaller while retaining the basic cell design illustrated 
in Fig. [T] there are two things that can be considered: 
increasing the cross sectional area of the transfer tube 
and decreasing the length of the transfer tube. We note 
that with the large-volume cells that are more resistant 
to beam current, design constraints make it difficult to 
significantly shorten the transfer tube, and they are cur- 
rently at most 50% longer than the shortest lengths used 
in early two-chambered targets such as those employed 
during E142. We will thus focus here on the cross sec- 
tional area. 

To better understand the limitations of single transfer- 
tube configurations for future targets, we consider a hy- 
pothetical target design for an approved experiment at 
JLab (E12-06-016) that will measure G^ up to Q^ = 
lOGeV^ [^. The experiment will run at 60 ^A with 
a target length of 60 cm instead of 40 cm, a luminos- 
ity equivalent to running 90 /iA into the target of Fig. [T] 
Here we will assume a single transfer tube is employed 
instead of the convection-based design that is actually 
specified in the proposal. 

The single best way to make a target less susceptible 
to beam current is to make the cell larger and use more 
lasers, thus making relaxation from the beam a smaller 
perturbation to the target as a whole. The proposal for 
E12-06-016 calls for a target containing roughly 7 STP 
liters of gas, in contrast to the roughly 3 STP liters con- 
tained in the targets illustrated by Fig. [l] Using Fig. [l] 
as a starting point, we will assume a transfer tube diam- 
eter equal to the diameter of the target chamber (about 



III. CONVECTION DRIVEN CELLS 

We describe next a variant of the target cell geome- 
try depicted in Figure 1. There are still two chambers, 
a pumping chamber and a target chamber, but the two 
chambers are connected by two transfer tubes instead of 
one. With this design, it is possible to induce convec- 
tion, thus causing rapid transfer of gas between the two 
chambers. Furthermore, all that is required to induce 
convection is to maintain a temperature differential be- 
tween the vertical segments of the two transfer tubes. By 
controlling the temperature differential, the speed of the 
convection can be adjusted. With rapid mixing of gas be- 
tween the two chambers, the aforementioned polarization 
gradients can be made negligible, even if the distance be- 
tween the pumping and target chambers is substantially 
increased. 



A. Experimental Setup 

To demonstrate the feasibility of convection-driven po- 
larized ^He target cells, we have constructed a prototype 
with the geometry and dimensions illustrated in Fig. [4] 
The cell was constructed entirely out of aluminosilicate 
glass (GE 180), and was sealed after being filled with 7.2 
amg of ^He and a 0.11 amg of N2. The pumping cham- 
ber also contained several tens of milligrams of a hybrid 
mixture of potassium and rubidium. 

The '^He gas in the new prototype was polarized in the 
same manner as in our other target cells. The pumping 
chamber was surrounded by a forced-hot-air oven, con- 
structed largely out of ceramic and glass. While polar- 
izing the "^He, the oven was maintained at temperatures 
that were typically between 200 — 235° C, resulting in a 
vapor pressure of alkali-metal atoms corresponding to a 



number density on the order of lO^^cm"^. The rubidium 
atoms were opticahy pumped using laser hght from high- 
power diode-laser arrays with a wavelength of 795 nm, 
and as described in Refs. [T3| and [T5], quickly shared 
their polarization with the potassium atoms that were 
also present. Subsequent spin-exchange collisions with 
the ^He atoms resulted in the buildup of substantial nu- 
clear polarization. 

The temperature differential used to induce convection 
was maintained by using a second forced-hot-air "con- 
vection heater" , installed on one of the transfer tubes as 
illustrated on Fig. [4] (we note that, in reality, the convec- 
tion heater also covered much of the horizontal portion of 
the transfer tube, but that this portion contributed neg- 
ligibly to the speed of the convection). With a portion 
of one of the transfer tubes at an elevated temperature, 
the gas contained therein had a lower density than the 
corresponding gas in the other transfer tube, and thus 
experienced a small buoyant force which induced convec- 
tion. By controlling the temperature of the convection 
heater, the gas flow could be controlled in a stable and 
reproducible fashion. 

The flow of the gas was monitored using an NMR tag- 
ging technique. A "slug" of gas within a small section 
of one of the transfer tubes was depolarized by subject- 
ing it to a pulse of RF tuned to the Larmor frequency of 
the "^He nuclei. The RF was delivered using a small coil 
(labeled in Fig. H as the "Zapper coil" ) wrapped directly 
around one of the transfer tubes. NMR signals were then 
detected at each of four locations along the target cham- 
ber as indicated on Fig. [4] using small "pickup coils" . 
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FIG. 4: Prototype convection-based target cell. The pumping 
chamber is placed inside of an optical pumping oven. The 
right transfer tube is heated while the left transfer tube is kept 
at room temperature. The two transfer tubes have different 
densities which creates a counter-clockwise convection current 
in the cell. The zapper coil is used to depolarize a slug of gas. 
This slug is then monitored as it travels through the pickup 
coils on the target chamber. We note that a small horizontal 
portion of the transfer tube was also heated, but for clarity is 
not shown, since it did not contribute to driving convection. 



The movement of the slug of gas could then be tracked 
by monitoring NMR signals from each of the four pickup 
coils. Signals from the pickup coils were obtained once 
every two seconds using the NMR technique of adiabatic 
fast passage (AFP) [21] . Representative examples of such 
signals are plotted in Fig. [5] as a function of time. Time 
zero in this plot corresponds to the moment when a slug 
of gas was tagged. 
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FIG. 5: Data used to visualize gas flow are shown in which 
the NMR signal from four pickup coils are plotted versus time. 
The oven temperature was 215°C, and the transfer tubes were 
24° C and 50° C respectively. The data indicate a gas flow 
velocity of 20 cm/min in the target chamber. 

It is readily apparent from Fig. [5] that a transient dip 
occurs in each of the signals from the four pickup coils. 
This dip corresponds to the passage of the depolarized 
slug of gas, and it can be seen that the transient occurs 
at successively later times for each of coils 1-4. Given 
the known positions of the regularly spaced pickup coils, 
the difference in time between the transients associated 
with each of the four pickup coils provides a measure of 
the speed with which the tagged slug of gas was moving. 
The measurement illustrated in Fig. [5] corresponds to a 
target-chamber gas velocity of 20 cm/minute. 

Also apparent from Fig. [5] is the fact that each succes- 
sive dip becomes wider and more shallow. This is due in 
part to the fact that the gas flow within the cell is char- 
acterized by the classic parabolic Hagen-Pouiselle veloc- 
ity distribution. That is, the velocity as a function of 
the distance r from the middle of the tube has the func- 
tional form v{r) = Wmax(l — r^/i?^). Hence, the slug of 
gas, which is initially fairly localized to the region around 
the "zapper coil", becomes increasingly spread out as it 
moves through the target chamber. Additional spread- 
ing also occurs because of diffusion. We note also that 
the NMR signal decreases as a function of time. This 
is due largely to polarization losses that occur with each 
AFP measurement. While the loss from each individ- 
ual measurement is quite small (on the order of 1%), the 
accumulation of many such losses is quite substantial. 



9 



We note that in the normal operation of a polarized ^He 
target measurements are typically made not every two 
seconds, but rather every two to four hours. 



B. Temperature Dependence of the Gas Velocity 

Hagen-Pouiselle flow occurs whenever a pressure dif- 
ferential between two ends of a pipe causes the laminar 
flow of a viscous gas or fluid [32 El]- In equilibrium, 
the driving force from the pressure differential i^driving 
must be equal to the retarding force -^retarding from the 
viscosity: 



-^driving — -^retarding 



(35) 



For the case of a pipe that is circular in cross section, 
Eqn. [35] must be satisfied for each annular ring of fluid 
of thickness dr, a condition which leads to the equation 

AP2nrdr = -2TTnl— ( r-^ ) dr , (36) 

dr \ dr J 

where AP is the pressure differential, rj is the viscosity 
of the fluid and I is the length of the pipe. Imposing 
the boundary condition that the velocity of flow must go 
to zero at the perimeter of the pipe, the solution to this 
differential equation is 



v{r) 



1 AP(r2 - i?2) 
4 ^l 



(37) 



where R is the radius of the pipe. It is the velocity dis- 
tribution given by Eqn. [37| that is often referred to as 
Hagen-Pouiselle Flow. 

In the case of our convection cell the driving force is 
due to the small buoyancy force that results from main- 
taining a vertical portion of one transfer tube at a higher 
temperature than the corresponding section of the other 
transfer tube: 



-Fv 



buoyancy 



^ApVtg 



(38) 



where Vt is the volume of the vertical portion of the trans- 
fer tube that is being heated, Ap is the difference between 
the average densities of the heated and unheated por- 
tions of the transfer tubes, and g is the acceleration due 
to gravity. We can express Ap as follows: 



Ap = pTc 



1 



(39) 



where p is the density of the gas in those portions of the 
cell that are not heated, Tq is the temperature of those 
portions of the cell that are not heated, and Th is the 
temperature of the portion of the transfer tube that is 
being heated (both temperatures are in Kelvin) . For the 
case being discussed here, the pressure differential that 
appears in the left-hand side of Eqn. [36] is thus given by 



AP = PbuoyancyMt = Ap/ig, 



(40) 



where At is the cross sectional area of the transfer tube, 
and h is the length of the portion of the transfer tube 
that is heated. 

If our convection cell could be treated as a long straight 
tube, we could simply substitute Eqns. [39] and [40] into 
Eqn. |37] to obtain an expression for the velocity. Our 
convection cells are more complex, however which com- 
plicates the expression that appears on the right-hand 
side of Eqn. 36 There are multiple sections of tubing, 
each with its own radius, as well as bends etc. The veloc- 
ity of the gas will be different in each section, and even 
the viscosity will be different depending on the tempera- 
ture. Luckily, however, as will be shown in the appendix, 
the continuity equation ensures that the velocity in each 
section is related in a simple linear fashion to the velocity 
in the other sections, so it is still possible to solve Eqn. [36] 
exactly. In essence, the quantity rjl that appears on the 
right-hand side of Eqn. [36] must be replaced by a single 
quantity fc^£ , which still has the dimensions of viscosity 
times length, but incorporates the full complexity of the 
cell. The solution can accordingly be written in the form 



(r^ - R^) 
v{r,Tii) = —- phgTc 



4fc, 



»,£ 



1 



1 



Tr. Th 



(41) 



where r is the radial coordinate in the target chamber, 
and R is the radius of the target chamber. The temper- 
ature dependence of v is largely dominated by the factor 
(t — 7~)- "^^^ quantity fc^£, however, is also dependent 
on temperature, although for the range of values of Tjj 
that we consider, the temperature dependence of k,^c on 
Th is relatively weak. 

It is important to understand the relationship between 
the observed velocity of the gas, v°^^, as indicated by data 
of the sort sh own in Fig. [5| and the velocity distribution 
given by Eqn. 41 The natural way to compute v°^^ is by 
taking the physical separation of adjacent pickup coils, 
and dividing by the separation in time between the min- 
ima of the corresponding transients. If we consider the 
limit in which the distance between the zapper coil and 
the pickup coils is long compared to the length of the 
zapper coil itself, it is straightforward to show that, to 
a good approximation, the above method of computing 
corresponds to the maximum value of the velocity 



,,obs 



given by Eqn. |41[ 



which results from setting r = 0. 



To a good approximation, we can express u" 
simplified form 



A 



l + PiAT 



1 
Tc 



1 



the 



(42) 



where all quantities not dependent on Th have been ab- 
sorbed into the constant A, and the temperature depen- 
dence of kjjc is accounted for by the factor 1 + /3iAT, 
where AT = Th — Tq. For the conditions we have con- 
sidered, /3i is on the order of 10~'^/° C. 

In Fig. [6]we have plotted v°^^ as a function of Th for a 
range of temperatures. For each point, the velocity was 
computed using data such as those shown in Fig. [5] We 
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also show in Fig. [6] with a sohd black line a fit of the 
data to a function of the form of Eqn. |42j The qual- 
ity of the fit is clearly quite good, and yields the val- 
ues A = 7 AT [22) X lO^K-cm/min, Tcom = 24.3(8), and 
/3i = —0.2(2) X 10"'^/° C. For comparison, using our 
best knowledge of the cell geometry and densities, with 
Tcoid = 24.5°C, we compute A = 9.14 x lO^K-cm/min 
and Pi = 1.03 x 10" V° C (see Appendix A). Given the 
uncertainties of some of the quantities with which we are 
working, particularly in describing the retarding forces 
associated with our relatively complicated cell geometry, 
this agreement is quite reasonable. More importantly, 
the agreement is more than sufficient to suggest that we 
have an acceptable quantitative understanding of the pa- 
rameters infiuencing the convective flow from a practical 
perspective. 




Heated Transfer Tube Temperature 



FIG. 6: The velocity of the gas in the target chamber of the 
convection-based cell is shown as a function of the temper- 
ature T of the heated transfer tube. The oven temperature 
was 215°C and the unheated transfer tube was at 24°C. 



C. Convection Transfer Rates and the Elimination 
of Polarization Gradients 

Ultimately, the value of a convection-driven target cell 
is measured by the degree to which polarization gradients 
can be avoided between the pumping chamber and the 
target chamber. It is critical that, as gas is depolarized 
by an electron beam, freshly polarized gas is delivered 
from the pumping chamber sufficiently quickly. Ideally, 
one would like the ratio of the polarizations of the two 
chambers to be as close to unity as possible. Fortunately, 
when convection is used to transfer gas between cham- 
bers, the polarization gradient is suppressed because the 
transfer rates are quite high. 

For the range of velocities measured in Fig. l6J we find 
values for dtc in the range of 4.9 - 81 hrs~^. Even with 
relatively fast relaxation in the target chamber (for exam- 
ple, consider Ftc = 1/10 hrs~^), the polarization gradient 



(Eqn. p}0[ ) will be very small in the presence of such fast 
transfer rates (A < 0.02, or Ptc/Ppc > 0.98). 

In Fig. [2J the ^He polarization in both chambers of a 
traditional target cell is plotted as a function of time. 
The ratio of target chamber to pumping chamber polar- 
ization for this plot begins at zero, and gradually climbs 
to a value significantly less than unity, reflecting a sub- 
stantial polarization gradient. In Fig.lTJ we plot this ratio 
for a convection-driven cell for three different operating 
conditions. In all cases, the temperature of the oven was 
held at 215 °C. The three curves correspond to different 
temperatures Th of the heated transfer tube. For the 
data shown with the open squares, the transfer-tube set 
temperature was 24 °C, the same as the other (unheated) 
transfer tube. This case corresponds to no driven con- 
vection, and is associated with a polarization gradient 
of 8%. For the data shown with the filled triangles and 
the filled circles, the set temperatures were 50 °C and 
100°C respectively. These two conditions corresponded 
to target-chamber gas velocities of approximately 19.9 
cm/min and 48.5 cm/min. In both of these cases, the 
ratio of the polarizations of the target chamber and the 
pumping chamber quickly reached a value approaching 
unity. It is notable that there is very little difference be- 
tween these last two curves despite substantially different 
gas velocities. In short, as soon as convection rather than 
diffusion is responsible for the gas transfer between the 
two chambers, polarization is relatively uniform through- 
out the cell. 
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FIG. 7: The ratio of the polarizations of the target chamber 
and the pumping chamber, Ptc/Ppc, is shown versus time for 
3 spin-ups corresponding to different convection velocities. 



IV. OUTLOOK FOR FUTURE TARGETS 

With advances in SEOF, it is now possible to run po- 
larized '^He targets at significantly higher luminosities. 
Using a single transfer tube, however, even currently ap- 
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proved experiments would suffer polarization gradients of 
12% or worse, and subsequent experiments at even higher 
luminosities would be even more limited. The work de- 
scribed here demonstrates that with only minor changes 
to current designs, convection-based cells can be built in 
which polarization gradients are extremely small. 

In many ways, however, the most compelling reason 
to adopt convection-based cells is the enormous flexibil- 
ity that is gained in being able to choose the distance 
between the pumping and target chambers. Experience 
has shown that the probability of a cell rupturing goes 
up markedly after 4-6 weeks exposure to something on 
the order of 12 /xA of beam. Most commonly, the part of 
the cell that breaks is the pumping chamber. Stress in a 
sphere goes up linearly with the radius. This makes the 
relatively large pumping chamber particularly vulnerable 
to radiation damage. With convection-based cells the 
distance between the pumping chamber and the target 
chamber can be greatly increased with minimal implica- 
tions for performance. Among other things, it becomes 
practical to incorporate radiation shielding for the pump- 
ing chamber. 

There are at least two additional reasons that flexibil- 
ity in the configuration of target cells is highly desirable. 
As discussed earlier, the primary way to make targets 
more tolerant of beam is to make them larger. The more 
gas that is being polarized per unit time, the less relevant 
it is that some small portion of that gas is being depolar- 
ized by the beam. It is not practical, however, to make 
the pumping chamber arbitrarily large. As the stress in 
the glass walls increases, the cell becomes more prone to 
rupturing, even if the wall thickness is increased. Fur- 
thermore, the intensity of the optical pumping radiation 
goes like total power over the square of the radius while 
the volume goes like the cube of the radius. Hence, since 
the required laser power scales with volume, the required 
intensity of light scales linearly with radius. There are 
already indications that we are near a threshold at which 
laser intensity damages the cells from the perspective of 
spin relaxation. If multiple pumping chambers were used 
instead of one, it would become possible to limit both 
laser intensity and the size of each individual pumping 
chamber. If diffusion were the only mechanism for mov- 
ing polarized gas inside the target, the implementation of 
multiple pumping chambers would be quite challenging. 

We close by citing a concrete example that underscores 
the desirability of adopting convection-based target cells. 
We show in Fig. [8] results from a bench test of the cell 
"Brady" (discussed earlier) in which the primary goal 
was to achieve the highest polarization possible. Brady, 
like the the targets used for the G^ measurements [3], 
contained an alkali-hybrid mixture of K and Rb. Unlike 
during the G^ experiment, however, the results shown 
in Fig. [8] were obtained using diode-laser arrays that 
were spectrally-narrowed such that their line widths were 
about ten times narrower than their broadband counter- 
parts. The polarization in this test saturated just over 
70%. Brady was one of two target cells used to collect 



data on single-spin asymmetries in semi-inclusive deep 
inelastic scattering, and was similar in design, although 
slightly smaller, than the configuration shown in Fig. IT] 
During the experiment, the in-beam polarization of the 
■^He averaged around 55% [35,, despite the fact that the 
cell-averaged polarization, averaged over the entire ex- 
periment, was over 60%. If the targets used for this ex- 
periment had had convection-based gas mixing (which 
was just being explored at the time), significant gains in 
performance would have resulted. 
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FIG. 8: Shown is polarization as a function of time for the 
cell "Brady" . A polarization of 70% was ultimately reached. 
The set of tests giving rise to these data established the high- 
est polarizations observed to date among large polarized ^He 
target cells intended for electron scattering. 

In conclusion, with the advances in target perfor- 
mance brought about by a combination of alkali-hybrid 
spin-exchange optical pumping and spectrally-narrowed 
diode-laser arrays, the point has been reached when it is 
desirable to explore target-cell geometries that go beyond 
that shown in Fig. 1. We suggest that the convection- 
based gas mixing demonstrated in this work can play 
an important role in implementing high-luminosity po- 
larized •^He targets in the future. Not only will the use 
of convection reduce polarization gradients, it will also 
make possible greater flexibility in the placement of the 
pumping chamber with respect to the target chamber. 
Among other things, this flexibility will enable the use of 
radiation shielding, and even the use of multiple pumping 
chambers. Convection-based polarized "^He target cells 
open new possibilities that collectively address several of 
the important challenges facing the next generation of 
high-luminosity polarized '^He targets. 
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Appendix A: Estimate of the magnitude of flow in 
the convection cells. 



Here we formulate estimates for the parameters that 
appear in Eqn. |42[ 



where h is the vertical length of transfer tube that is held 
at an elevated temperature. In this model we approxi- 
mate the pumping chamber as a cylinder with transfer 
tubes entering axially, and identify five distinct regions in 
the cell as is indicated in Fig. [9] Each region is identified 
as a pipe of length Zj, radius i?i, cross sectional area Ai 
and temperature Tj. We further assume that Ti — T^, 
Tg = Tg and i?i = i?2 = i?3- Finally, we identify Tc = Ti 
and Th = T2. We note that both the density and viscos- 
ity of the gas are temperature dependent. 

The continuity equation, pjAjVj = piAiVi and some 
distance rescaling provide further simplification: 



1. The Viscosity of ^He 

Since our temperatures are in the classical regime 
(T ^ 3 K), the viscosity ryncs of '^He can be calculated 
from the viscosity ?7hc4 of *He using [36 : 



f7Ho3 



™Hc3 



-?7Ho4 
mHo4 

0.8681 ma 



(Al) 
(A2) 



We parameterize the ^He viscosity in the range of 0—300° 
C following Kestin et al. |37) . 



r7Ho4 = A + BxT + CxT^ 

where T has units of ° C and 

A = 18.82(2) /iPa-s , 
B = 0.0456(2) ^Pa-s/°C 
and C = -13.8(6)pPa-s/(°C)2 



(A3) 



(A4) 
(A5) 
(A6) 



At 20°C, r^Hcs = 17.12^Pa-s. 

The flow in a pipe is Laminar if the Reynold's number 
is below 2300 [38]. The Reynolds number is defined as 



i?e = 



2Rpv 



(A7) 



where p is the density of the fluid. A pipe 1.2cm wide 
that is filled with 8 amagats '^He will have laminar flow 
at 20° C (p w 1kg /m^) provided v < 20000 cm/min. 



2. Flow in the Convection Cell 

The flow in the convection cell arises from a forced 
density difference between the two transfer tubes - one 
tube is maintained at room temperature while the other 
is heated (see Fig.ji]). We modeled the convection cell as 
five contiguous pipes as is illustrated in Fig. [9] Eqn. 36 
becomes 
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Apgh2Trrdr = —2Try Vih^r' 



^^Sl'^^' 



(A8) 



Vi 



and 



dr-i 



PiRl 
PiRf 

Ri 

— fi 
Ri 

Ri d 
Ri dri 



Vi 



(A9) 
(AlO) 
(All) 



Since vi and r^ have been expressed in terms of vi and 



■^^^^^^^^^^ 







^^ 



FIG. 9: Diagram indicating the different temperature regions 
used to describe the convection cell. 

ri, we'll drop their subscripts. Finally, 



v{r) 



Apgh I 



Ri) 



Vi Ml +'4-52- +mh 



PiRi 
■7^ 



m 



53 V 



+ / 



where v{r) is the velocity in region 1. 



(A12) 



Eqn. (A12) assumes that there are no "minor losses" 



in the system, where a minor loss represents a pressure 
drop due to a sudden change in flow from a pipe fitting 
or a pipe expansion or contraction; the term "minor loss" 
refers to a loss that is small relative to the overall length 
of pipe under consideration 33, 39 . In our case, due to 
the relatively short length of the cell, the minor losses 
are actually significant. The actual cell, as illustrated in 
Fig.|4j consists of four elbow bends, two tees, and four ex- 
pansions/contractions. The retarding forces these losses 
exert on the gas can be approximated by considering in- 
stead an equivalent length of straight pipe. 

The expansions/contractions have a relatively small 
loss, which is equivalent in magnitude to the loss that 



13 



would be incurred passing through a pipe of length 

2RK 



L 



equivalent 



/ 



(A13) 



where R is the tube radius, / is the friction factor 
(/ = 64/i?e for laminar flow) and K is the fluid- 
independant resistance coefficient. For sudden expan- 
sions/contractions, 



K 



expansion 



and K, 



contraction 




(AM) 



(A15) 



Gas flowing through the transfer tube/target chamber 



0.18, K, 



expansion 



0.35. 



junction has it'contraction 

At u = 60cm/min {Re « 10), this gives a negligible 
^equivalent ~ O.lcm - gas flowing between the pumping 
chamber and the transfer tube will have an even smaller 

-'^equivalent ■ 

The losses in the bends, however, are much greater. 
We model the loss coefficient in the bends using the 3-K 
method of Darby [ID] , 



K 



Re 



+ kAi 



£,0.3 



(A16) 



where D is the diameter of the pipe (in inches) and 
Ki,Ki,Kd are geometry-dependent loss coefficients. We 
approximate our glass bends as flanged, welded bends 
with rb/D = 2 (here, rt, is the radius of the bend); such 
bends have Ki = 800, Ki = 0.056, Kj, = 3.9. For laminar 
flow, the 3-K method gives 



-'^equivalent 



2R 
64 



Ki + ReKi 1 



Kd 

£)0.3 



(A17) 



We treat the transfer tube/target chamber tee junctions 
as elbows (effectively ignoring the dead-end branch of the 
tee) . The system therefore has five bends in temperature- 
region 1 (which have a total equivalent length of approx- 
imately 63 cm) and one bend in temperature-region 2 
(which has an equivalent length of approximately 13 cm). 
Finally, using Eqn. |39|for Ap we find 



v(r) 



where 



l{r^~Rl)phgTc{^ 



1 



vi[i[ + k'-^)+mi'.^^ 



m^ 



H' 



+ 1 



I'l 



li + 5x 



2Ri 



Ki+ReiKi 1 



Kd 

^0.3 



I2 + 



2R2 

64 



Ki + Re^Ki 1 



Kd 
r,o.3 



(A18) 

(A19) 
(A20) 



Note that the Reynold's number is dependent on the ve- 
locity of the gas. 



Table \K2\ lists values for R and I. Measurements of R 
(which is the inner diameter of the tube) require a knowl- 
edge of the thickness of the glass tube. We measured the 
thickness of the glass by observing interference patterns 
using a scannable single-frequency laser. Using this in- 
formation, Eqn. I Al8| predicts velocities that agree within 
20% of the measured value (see Fig. l6| . 



Ri 


R2 


R3 


R4 


R5 


h 


li 


I2 


I3 


I4 


I5 


.498 


.521 


.502 


.806 


4.034 


4.76 


24.99 


20.87 


15.18 


40.32 


5.6 



TABLE I: Best Guesses for Convection Cell Dimensions (cm) 



The maximum velocity in Eqn. AI8 (corresponding to 
r = 0) can be written as 



V Tcoldtt Thottt 7 

1 + (3{AT) 



where 



AT 
A' 

P 

kriC 



RJphgTc 
^kr,c ' 

/3iAr + /32(AT)2 



V2I 



kiATf 

,piRl 



+V3 



T3 



P2RI 

r>2 d2 



AT = T2- Ti 

Pi 



,Rl,m 



■■Ri^T, 



7-.2 

0.868n; ^ 



and 
/33 = 0.8681/: 



i?2 



0.8681(5 

2C(ri-273))] , 
B + 2C(Ti - 273) 



C 



I Rj 
Rl 



c 



(A21) 

(A22) 
(A23) 

(A24) 

,(A25) 
(A26) 

(A27) 
(A28) 
(A29) 
(A30) 



In the above equations, all temperatures are in Kelvin. 
To describe the velocity in region 4 (the target chamber) 
instead of in region 1, we use Eqn. A9 and replace A' 

with A = i-^)A' . This is the quantity that appears 

"4 

in Eqn. |42] that characterizes the magnitude of the gas 

velocities shown in Fig. [6] 

Evaluating the above equations in terms of our best 

guesses for cell dimensions and temperatures, with Tq — 

24.5°C gives A = 9.14 x lO-^K-cm/min, /3i = 1.03 x 10^^^ 

P2 = 1.26 X 10^6, and h = -4-26 x 10"^°. It is clear 

from these values that it is quite reasonable to neglect 

the terms involving /32 and j3^, which results in Eqn. [42] 

from section III-B. 
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